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ABSTRACT 

Characterizing rotation, infall and accretion disks around high- mass protostars is an important topic in 
massive star formation research. With the Australia Telescope Compact Array and the Very Large Array 
we studied a massive disk candidate at high angular resolution in ammonia (NH3(4,4) & (5,5)) tracing 
the warm disk but not the envelope. The observations resolved at ^0.4" resolution (corresponding to 
^1400 AU) a velocity gradient indicative of rotation perpendicular to the molecular outflow. Assuming 
a Keplerian accretion disk, the estimated protostar-disk mass would be high, similar to the protostellar 
mass. Furthermore, the position- velocity diagram exhibits additional deviation from a Keplerian rotation 
profile which may be caused by infalling gas and/or a self-gravitating disk. Moreover, a large fraction 
of the rotating gas is at temperatures >100K, markedly different to typical low-mass accretion disks. 
In addition, we resolve a central double-lobe cm continuum structure perpendicular to the rotation. We 
identify this with an ionized, optically thick jet. 

Subject headings: stars: formation - stars: individual (IRDC 18089-1732) - stars: early-type - stars: 
rotation - stars: winds, outflows 



1. INTRODUCTION 

Outflow and jet studies in high-mass star-forming re- 
gions have revealed ample indirect evidence that accre- 
tion disks should be common in these regions. However 
direct observational c haracterization of such disks has so 
far b een difficult fe.g. JBeuther et aT]|2007l : [Cesaroni et alj 
l2007lh While this is partly caused by the typically large 
distances and the clustered mode of massive star forma- 
tion, an additional problem arises because these regions 
are still deeply embedded within their natal cores, they 
have visual extinctions of the order 1000, making it diffi- 
cult to disentangle the dense core emission from the gen- 
uine accretion disks. In the past, different molecular line 
tracers were employed for differe nt massive star-forming 
regions (e.g.. ICesaroni et aT1l2007t ). however, it is interest- 
ing to note that for many sources only one or two molecular 
lines exclusively indicated rotating motions whereas many 
other promising sp ectral lines showed no such signatures 
(e.g. . iBeutherl [20071 ). The problem of identifying a reliable 
massive disk tracer prohibited statistical studies of larger 
disk samples, until today. Here we present new observa- 
tions toward the massive disk candidate IRAS 18089-1732 
in the high-excitation ammonia lines NH3(4,4) and (5,5), 
with lower energy levels of 200 and 295 K, respectively. 
While the high excitation temperatures ensure that we are 
not tracing t he cold gas envelope b ut only the warm central 
regions (e.g.. ICesaroni et aT1ll998f ). recent radiative trans- 
fer calculations of 3D hydrodynamic massive core simu- 
lations indicated that the wavelength regime < 100 GHz 
should be pa rticularly advantageous to studying the inner 
disk regions (|Krumholz et al.ll2007af ). This work predicts 
high dust optical depths at frequencies > 100 GHz poten- 
tially limiting spectral line studies of the central disk in 
the (sub)mm regime. 

The target IRAS 18089-1732 has previously been stud- 
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ied at (sub)mm wavelengths with the Submillimeter Ar- 
ray, and a velocity gradient was identified from east to 
west, perpendicular to the molecular outflow emanat- 
ing f rom the core approximate ly in north-south direc- 
tion (jBeuther et alj 12004 I2005D . However, these data 
did not resolve the velocity pattern of this structure 
and hence did not allow to study its kinematic proper- 
ties. The region is located at ^3.6 kpc, it has a bolo- 
metric luminosity of about 10 L©, and it exhibits sev- 
eral other signs of early massive star f o rmation (e.g., 
I Walsh et all 119981: ISridharan et all 120021: iBeuther et alj 
2002 allbl: IWilliams et all 12004 120051 : iFuller et all 12005 ; 
lEdris et al.ll2007[ ). 



2. OBSERVATIONS 

2.1. The Autralia Telescope Compact Array (ATCA) 

IRAS 18089-1732 was observed in April 2007 with the 
ATCA in a 1.5 km baseline configuration including an- 
tenna 6. The phase reference center was R.A. (J2000) 
18Hl m 51 s .4, Decl. (J2000) -17°31'28".5. We observed 
the NHa(4,4) and (5,5) inversion lines with the frequen- 
cies of the main hyperfine components at 24.139 and 
24.533 GHz, respectively. The velocity relative to the lo- 
cal standard of rest (i>i sr ) is ~ 33.8 km s -1 . Phase and 
amplitude were calibrated by regular observations of the 
quasar 1829-207. Bandpass and flux were calibrated with 
observations of 1741-038 and 1934-638. The spectral res- 
olution of the observations was 62 kHz corresponding to a 
velocity resolution of ~ 0.8 km s -1 . Applying a robust 
weighting of 0.5, the la rms per 0.8kms _1 channel is 
~2.5mJybeam _1 . The synthesized beam is 2.8" x 0.4" 
with a position angle of 8.5°. Because the o utflow is in 
the north-south direction (jBeuther et al.ll2b04h and the ex- 
pected accretion disk in the east-west direction, the small 
synthesized beam along the east-west axis is fortunate for 
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the analysis. 

2.2. The Very Large Array (VLA) 

To improve the beam shape, we observed the same re- 
gion in October 2007 with the VLA/eVLA hybrid array 
in an extended configurations (BnA) in the NH 3 (5,5) line. 
Due to mediocre weather conditions and a worse perfor- 
mance of the retrofitted eVLA antennas at high frequen- 
cies, we could not usefully include the longest baselines, 
restricting the data to uv-ranges <500kA. Bandpass and 
flux calibration were conducted with 1331+305 and 1733- 
130. Fast switching between source and phase calibrator 
1832-105 were applied to calibrate the phases and ampli- 
tudes. The original spectral resolution of 48.8 kHz was 
smoothed in the final data-cube to l.Okms -1 . Applying 
natural weighting and restoring the data with a circular 
beam (degraded to the larger axis of the fit to the dirty 
beam) the synthesized beam of these data is 0.47". It 
was also possible to extract a 1.2 cm continuum map from 
the line-free part of the spectrum. For this image we in- 
cluded longer baselines as well and the synthesized beam 
is 0.26" x 0.19" with a PA of 84° from north. 

3. RESULTS AND DISCUSSION 

Figure [T] presents a compilation of the ATCA data com- 
pared to the submm c ontinuum emission tr acing the cen- 
tral dust and gas core (jBeuther et al.|[2uu5t ). As expected 
the warm gas observed in the NH3 lines traces the cen- 
tral core well. However, more interesting is the veloc- 
ity structure of the region. The two first moment maps 
(the intensity weighted velocities) of the NHa(4,4) and 
(5,5) data-cubes both exhibit a velocity gradient from east 
to west, perpendicular to the molecular o utflow emanat- 
ing a pproximately in north-south direction dBeuther et alj 
2004). The north-south elongation seen in Fig. Q] is caused 
by the non-symmetric beam shape. To improve on that, 
we re-observed the source in the NH3(5,5) line with the 
VLA, Fig. [5] presents and integrated and a first moment 
map of these data. The velocity gradient in east-west di- 
rection is obvious again. The fact that the NH3 structure 
is not flattened could be either due to a viewing-angle of 
the rotating structure close to face-on, or, since the NH3 
features in the north and south are close to the V\ sr , they 
may also contain contributions from the ambient core. 

In addition to the spectral line data, Fig. [5] presents the 
simultaneously with the VLA observed 1.2 cm continuum 
map. The cm emission shows a double- lobe structure in 
north-south direction with the dip centered on the submm 
continuum peak. The integrated emission of this double- 
lobe is 7.4 mjy In combination with the previously ob- 
served 3.6 cm flux of 0.9 mjy, the spectral index between 
these two waveleng ths is ~2, cons istent with an ionized, 
optically thick jet (|Revnoldsl [i~986) . Since the molecular 
outflow is also in north-south direction we identify this 
double-lobe cm continuum emission with the central jet 
perpendicular to the rotation axis of the gas. 

To get a better impression of the velocity field along 
the inferred rotation axis, Fig. [3] shows a position- velocity 
digram cut through the submm continuum peak posi- 
tion along the east-west axis. For this analysis, we use 
the ATCA NH3(4,4) data because they combine the best 
signal-to-noise ratio with the highest spectral and angular 



resolution along the given axis. The velocity gradient is 
clearly depicted here as well. With a spatial resolution of 
0.4" along this axis, the structure with an approximate size 
of 2" is well resolved and consistent with rotation of a po- 
tential massive accretion disk. Furthermore, with a linear 
resolution of ~1400AU, these observations allow us, for 
the first time, to better characterize the kinematics of the 
rotating structure. 

Refering to the relatively simple picture of low-mass 
accretion disks wh i ch are in Kepleria n rotation (e.g., 
lOhashi et all Il997t iDutrev et all l2007t ). we can calcu- 
late the rotationally supported binding mass M rot for 
IRAS 18089-1732 assuming equilibrium between the rota- 
tional and gravitational force at the outer radius of the 
disk. 



=> M rot [M ] = 1.13 10~ 3 x Sv 2 [km/s\ x r[AU] (2) 
Here, r is the disk radius (~ 1" ~3600AU), and Sv 
half the velocity regime observed in the first moment maps 
(~3kms _1 , Fig. [J). Equations[T]&[2]have to be divided by 
sin 2 (i) where i is the unknown inclination angle between 
the disk plane and the plane of the sky (i = 90° for an 
edge-on system). With the given values, we can estimate 
M ro t to ~ 37/(sin 2 (i)) M Q . This mass estimate is com- 
parable to the core mass estimates from the submm dust 
continuum emission (about 45 M ? . iBeuther et al.l 12005'). 

A different way to assess potential (non-)Keplerian com- 
ponents in the velocity structure is to draw an expected 
Keplerian velocity profile for a 15 M protostar on the 
position- velocity diagram (Fig. [31 for the mass estimate see 
the following paragraph) . While for blue-shifted velocities 
and positive offsets the position- velocity structure does ap- 
proximately follow the Keplerian profile, the red-shifted 
emission at negative velocities deviates more strongly from 
the synthetic profile, even showing super-Keplerian veloc- 
ities in excess of the Keplerian speed. Furthermore, there 
are additional components in the position- velocity digram, 
for example the strong peak at (31kms _1 , —0.2") that 
does not fit a typical Keplerian disk. How can we explain 
the non-Keplerian signatures in the rotating structure? 

Even at the highest angular resolution (continuum and 
line data in Fig. [T]), the source remains a single-peaked 
structure without exhibiting signs of multiplicity. Al- 
though massive star formation usually proceeds in a clus- 
tered mode, this indicates that most of the bolometric lu- 
minosity of 10 4 ' 5 Lq likely stems from the most massive 
central object. Assuming a main sequence star this would 
give an upper mass limit of 20 M Q for the central object. 
Since the source is likely still in its accretion phase, it 
has additional accretio n luminosity (jMcKee fc Tan| [2003; 
iKrumholz et al.|[2007b( ) and the actual protostellar mass 
may be lower. However, estimating for example the ac- 
cretion luminosity L acc = GM*M acc / for a star of 
M* = 15 Mq with a radius of i?» = 10 R© and an accretion 
rate of M acc = 10 -4 M Q yr _1 , one derives an L acc of only 
~5xl0 3 L Q . Although especiall y M a , rr and R* could both 
vary by a factor of a few (e.g., IKrumholz "eTaIll2007bf ) it 
is unlikely that the accretion luminosity contributes more 
than a factor two to the overall luminosity. Therefore, 
we estimate the mass of the central protostar to be likely 
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between 15 and 20 M Q . Compared to the mass of the po- 
tential protostar/disk system calculated either from the 
dust continuum emission or estimated via the rotational 
support assumption, the actual protostellar mass is more 
than a factor 2 lower. Hence, a Keplerian supported disk, 
where the disk mass is negligible compared to the proto- 
stellar mass, is not feasible in such a system. Because of 
the velocity gradient perpendicular to the outflow (Figs. Q] 
& \2§ and the comparatively low velocities compared to 
molecular outflows, it is unlikely that these features are 
due to unbound motions, e.g., from the molecular outflow. 
The more likely interpretation is that we are dealing with 
a massive rotating structure not in Keplerian rotation. 

Candidates for the non-Keplerian contributions are a 
rotating and infalling envelope and a potentially self- 
gravitating disk structure. For example, an infalling and 
rotating envelope produces similar features in the position- 
velocity diagram like the blue-shifted peak at (Slkms -1 , 
-0.2") in Fig. 0] (see, e.g., Fig. 10 in lOhashi etallll997|) . 
Such rotating, infalling envelope structures could resem- 
ble the proposed larger-scale toroids that are observed 
toward a small number of even mo re luminous sources 
(|Cesaronill2005t ICesaroni et al.ll2007D . It is expected that 
such toroidal non-equilibrium structures should be grav- 
itationally unstable, maybe forming companion sources 
(|Kratter fc Matznej 1200a IKrumholz et alil2007bf) . which 
then could produ ce additional k i nemat ic sig nature s . In a 
similar direction, iKeto fc Wool ([20061 ) and iKetol (|2007t) 
modeled the accretion flow around high-mass stars: In 
their case, the gas orb its on ballist i c trajectories aroun d 
a point mass (see also lUlrichl Il976l : iTerebev et al.lll984[) . 
and the accretion flow appears quasi-spherical at large 
distances from the star, flattening to a disk at smaller 
radii due to conservation of angular momentum (see Fig. 
9 in IKeto fc W ood 2006) . Also noteworthy are the red- 
shifted super-Keplerian gas velocities exceeding the Ke- 
plerian speed. While sub-Keplerian disks are predicted 
for disks with significant additional sup port against gray - 
ity than rotation (e.g., magnetic field, iGalli et ail [20061 . 
super-Keplerian velocities are expected if the disk con- 
tributes a significant fraction of the mass to the protostar- 
disk system. In this case, the outer disk rotates faster 
because it feels the gravit ational force of t he pro tostar 
and the inner disk (e.g., IKrumholz et~ai1 (|2007bf ) and 
Krumholz priv. comm.). 

In addition to the kinematic analysis, ammonia 
is a well known thermom eter for molecular gas 
(Wal mslev fc Ungerechtsl Il983). However, because even 
these high-excitation lines still have high optical depths 
(r > 1) we cannot accurately fit the hyperfine struc- 
ture of the lines. Therefore, reasonable rotational tem- 
perature estimates are not feasible. Nevertheless, the 
detection of spectral lines with excitation temperatures 
> 200 K implies that the average gas temperatures should 
exceed 100 K, consistent with previous temperature esti- 
mates based on HCOOCH3 and CH3CN observations at 



(sub)mm wavelengths (<Beuther et al.ll2004 [2005h . While 
the previous (sub)mm temperature estimates were based 
on data that did not resolve the kinematics of the rotat- 
ing structure and therefore could be attributed to the disk 
and/or envelope, the NH 3 data clearly show that the cen- 
tral rotating structure harbors large amounts of warm gas 
at temperatures >100K. This is markedly different to typ- 
ical low-mass disks where the bulk of the gas is at low 
temperatures <30K (|Pietu et al.ll2007f) . 

4. CONCLUSIONS AND SUMMARY 

From an observational point of view, this work shows 
that high-excitation lines at 1.2 cm wavelengths are well 
suited to study the kinematics around massive proto- 
stars. While massive accretion disk research in recent 
years has mostly focused on (sub)mm interferometer ob- 
servations, this may not be the ideal wavelength regime 
because of high optica l depth in the innermost regions 
(jKrumholz et al.ll2007af) . Although there exist examples of 
NH 3 disk studies, these were usually d one in the low exci- 
tation lines NH 3 (1 1) and (2,2) (e.g.. IZhang et all 119981: 
iBeltran et alj [20061 ). Such low-energy line observations 
suffer from contributions of the cold envelope material, 
and high-excitation NH3 observations like those presented 
here may turn out to be one of the best tracers of the 
central accretion disks around high-mass protostars. Sta- 
tistically larger samples have to be investigated this way 
to base our knowledge on more solid ground. 

In summary, these observations disentangle the kine- 
matics, the mass and the temperature of the rotating 
structure around a high-mass protostar, identifying a hot 
massive accretion disk candidate. The kinematic signa- 
tures are different from typical Keplerian velocity profiles 
known from low-mass star formation. These discrepancies 
could be produced by larger-scale rotating and infalling 
toroids or envelopes as well as by self-gravitating struc- 
tures within the massive accretion disk. The estimated 
gas temperature is >100 K which is also significantly larger 
than average temperatures from low-mass accretion disks. 
While the angular resolution of 0.4" is very high, corre- 
sponding to a projected linear resolution of ^1400 AU, 
we are still not resolving the inner accretion disk regions. 
Notwithstanding the fact that the larger-scale disk-like 
structure observed here deviates from typical Keplerian 
disk structures, it may be possible that such structures 
exist on even smaller scales closer to the central protostar. 
Future observational tests include high-spatial-resolution 
observations with ALMA and the completed eVLA, as well 
as spatially and spectrally resolved observations of vibra- 
tionally excited CO lines observable at mid-infrared wave- 
lengths (e.g., CRIRES on VLT). 

We like to thank Fabian Walter and Hendrik Linz for 
their help with the VLA data. H.B. acknowledges financial 
support by the Emmy-Noether-Programm of the Deutsche 
Forschungsgemeinschaft (DFG, grant BE2578). 
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Fig. 1. — The two left panels present in solid con tours the NHa(4,4) a nd (5,5) emission integrated from 31 to 37kms — 1 . The grey-scale 
shows the corresponding 860 /im continuum emission Beuther et al. (2005). The continuum emission is contoured from 10 to 90% (step 10%) 
of the peak emission of 1.4 Jy beam' 1 . The NH 3 maps are contoured in 3<r steps with 3cr values of 3.75 and 3.3 mjy beam for the (4,4) and 
(5,5) integrated maps, respectively. The two right panels show the corresponding 1st moment maps of the NH3(4,4) and (5,5) spectral line 
data-cubes. The asterisks mark the peak position of the submm continuum emission, and the synthesized bea m is shown at the bottom-left 
of each panel. The full line outlines the direction of the outflow as measured in SiO(8-7) (Bcuthcr et al. 2005). 
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Fig. 2. — The left panel shows the VLA NH3(5,5) emission integrated from 31 to 37kms —1 and contoured from 10 to 90% (step 20%) of 
the peak emission. The right panel presents the corresponding 1st moment map contoured from 31.5 to SG.Skms" 1 (step lkms -1 ). The 
white- black dashed contours show the 1.2 cm continuum emissi on. Contour levels s tart at 4<r and continue in la steps of 0.28 mjy beam - 1 . 
The asterisks mark the position of the submm continuum peak (Bcuthcr ct al. 2005), and the synthesized beams are shown at the bottom-left 
(grey NH3 and dashed 1.2 cm emission). 









63 

77 


1 , 



30 35 40 

Velocity [km/s] 



Fig. 3. — Position-velocity diagram of the NHa(4,4) data through the submm peak position in east- west direction. The contouring starts 
at the 2<r level and continues in la' steps (la ~ 2.5 mjy beam -1 ). The dotted lines mark the central position and the systemic velocity. The 
full line shows a Keplerian rotation curve with a central mass of 15 Mq, which obviously fail to fit the data well. 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0712.0579vl 



